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We systematically analyze the correlations between the various leptonic and hadronic flavor 
^ , violating processes arising in SUSY Grand Unified Theories. Using the GUT-symmetric 

relations between the soft SUSY breaking parameters, we assess the impact of hadronic and 



^ " leptonic flavor observables on the SUSY sources of flavor violation. 



I. INTRODUCTION 

Supersymmetry (SUSY) Breaking (SB) remains one of the biggest issues in physics beyond the 
Standard Model (SM). In spite of various proposals 1], we still miss a realistic and theoretically 
satisfactory model of SB. 

Flavor violating processes have been instrumental in guiding us towards consistent SB models. 
Indeed, even in the absence of a well-defined SB mechanism and, hence, without a precise knowledge 
of the SUSY Lagrangian at the electroweak scale, it is still possible to make use of the Flavour 
Changing Neutral Current (FCNC) bounds to infer relevant constraints on the part of the SUSY 
soft breaking sector related to the sfermion mass matrices [2] . 

The model-independent method which is adopted is the so-called Mass-Insertion Approximation 



* Unite mixte du CNRS et de I'EP, UMR 7644. 



(MIA) Q]. In this approach, the experimental hmits lead to upper bounds on the parameters (or 
combinations of) 5{- = A{ /m|, where A{- is the flavor-violating off-diagonal entry appearing in 
the / = (u, d, I) sfermion mass matrices and m?j is the average sfermion mass. The mass insertions 
include the LL/LR/RL/RR types, according to the chirality of the corresponding SM fermions. 
Detailed bounds on the individual 5s have been derived by considering limits from various FCNC 
processes j^-{lo|. 

As long as one remains within the simple picture of the Minimal Supersymmetric Standard 
Model (MSSM), where quarks and leptons are unrelated, the hadronic and leptonic FCNC processes 
yield separate bounds on the corresponding d'^'s and (5^'s. 

The situation changes when one embeds the MSSM within a Grand Unified Theory (GUT). In 
a SUSY GUT, quarks and leptons sit in the same multiplets and are transformed into each other 
through GUT symmetry transformations. In supergravity theories, where the supersymmetry 
breaking (SB) mediation to the visible sector is gravitational, the structure of the soft terms is 
essentially dictated by the Kahler potential. If the effective supergravity Lagrangian is defined at a 
scale higher than the Grand Unification scale, the matter fields present in the Kahler function have 
to respect the underlying gauge symmetry which is the GUT group itself. Subsequently, the SB 
mediation would give rise to the usual soft terms which however now follow the GUT symmetry. 
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1^. 



Hence, we expect quark-lepton correlations among entries of the sfermion mass matrices 
In other words, the quark-lepton unification seeps also into the SUSY breaking soft sector. 

Imposition of a GUT symmetry on the soft SUSY breaking Lagrangian >Csoft entails relevant 
implications at the weak scale. This is because the flavor violating (FV) mass-insertions do not get 
strongly renormalized through Renormalization Group (RG) scaling from the GUT scale to the 
weak scale in the absence of new sources of flavor violation. On the other hand, if such new sources 
are present, for instance due to the presence of new neutrino Yukawa couplings in SUSY GUTs 
with a seesaw mechanism for neutrino masses, then one can compute the RG- induced effects in 

n 

terms of these new parameters. As has been noted earlier [lll j. even in such cases, the correlations 
between hadronic and leptonic flavor violating Mis survive at the weak scale as a function of these 
parameters. As for the flavor conserving (FC) mass insertions (i.e., the diagonal entries of the 
sfermion mass matrices), they get strongly renormalized but in a way which is RG computable. 

In conclusion, in SUSY GUTs where the soft SUSY breaking terms respect boundary condi- 
tions which are subject to the GUT symmetry to start with, we generally expect the presence of 
relations among the (bilinear and trilinear) scalar terms in the hadronic and leptonic sectors. Such 
relations hold true at the (superlarge) energy scale where the correct symmetry of the theory is 
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the GUT symmetry. After its breaking, the mentioned relations will undergo corrections which 
are computable through the appropriate RGE's which are related to the specific structure of the 
theory between the GUT and the electroweak scale (for instance, new Yukawa couplings due to the 
presence of right-handed (RH) neutrinos acting down to the RH neutrino mass scale, presence of 
a symmetry breaking chain with the appearance of new symmetries at intermediate scales, etc.). 
As a result of such a computable running, we can infer the correlations between the softly SUSY 
breaking hadronic and leptonic 5 terms at the low scale where we perform our FCNC tests. Ex- 
plicit examples of such correlations in the context of an SU(5) SUSY GUT will be provided in next 
Section. 

Given that a common SUSY soft breaking scalar term of /2soft at scales close to Mpianck can give 
rise to RG-induced 5'^'s and (5''s at the weak scale, one may envisage the possibility to make use of 
the FCNC constraints on such low-energy (5's to infer bounds on the soft breaking parameters of 
the original supergravity Lagrangian (£sugra)- Indeed, for each scalar soft parameter of i2sugra one 
can ascertain whether the hadronic or the corresponding leptonic bound at the weak scale yields 
the stronger constraint at the large scale. This constitutes the major goal of this work: we intend 
to go through an exhaustive list of the low-energy constraints on the various (J'^'s and 5''s and, 
then, after RG evolving such 5^s up to Mpiancki we will establish for each 6 of >Csugra which one 
between the hadronic and leptonic constraints is going to win, namely which provides the strongest 
constraint on the corresponding f^sugra* 

We will show that there exists a very interesting complementarity in the sensitivity of the soft 
breaking sector of >Csugra to the FCNC constraints provided by hadronic and leptonic physics. 

The second, related purpose of this paper is to fully exploit the common origin from single ^sugra's 
of hadronic and leptonic weak scale (5's to make use of leptonic (hadronic) FCNC constraints to 
limit 5'^'s (5''s). In other words, for instance in a SUSY SU(5) context, one can use a leptonic 
FCNC process like r — > /U7 to impose constraints on the {^Sq,^ j^-^ term of the hadronic sector which 
are stronger than those which are derived from genuine hadronic FCNC processes like 6 — > 57 

In this respect, the present work intends to provide a particularly striking example of the 
correlation between "low-energy" (i.e. weak scale) experiments and "high-energy" (i.e. GUT or 
Planck scale) SUSY Lagrangian. The effort in the coming years, provided LHC yields some SUSY 
evidence, will be to "reconstruct" the original supergravity Lagrangian from which our weak scale 
testable SUSY descends. As we know, such work of reconstruction will be rather painful if we 
have to rely only on LHC results; on the other hand, accompanying the SUSY direct searches at 
LHC with the powerful FCNC tests will prove to be quite efficient in our effort of tracing back the 
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original parameters entering the underlying supergravity Lagrangian. In the case of SUSY GUTs 
such role of FCNC processes is further enhanced because of the above mentioned hadron-lepton 
correlations. 

The paper is organized as follows. In the next Section we will provide an example of relations 
between hadronic and leptonic FC 6's at the weak scale within SUSY SU(5). In Section IIIH 
we discuss the impact of GUT breaking effects (both at the tree level and at 1-loop) and the 
subsequent basis dependence which seeps in to the GUT-symmetric relations. In Section IIVI we 
discuss our procedure for obtaining bounds detailing various constraints we have imposed on our 
supersymmetric spectrum. In Sections |V] and Wl\ we proceed with a separate analysis of the 
hadronic and leptonic 6's, to conclude with a study of their correlations in Section [VIII In Section 
IVIIH the conclusions include some kind of final "score" on the tightness of the bounds in GUT 
related 6's coming from the various FCNC processes and the outlook for the interplay of the kind 
of flavor physics we consider here with the possible LHC outcome. 

II. ON HADRON LEPTON FCNC RELATIONS IN SUSY GUTS 

In this Section we provide an example of the correlations between hadronic and leptonic 6's 
entering the weak scale MSSM Lagrangian when the underlying theory at the large scale, where 
supergravity lives, is restricted by a grand unified symmetry. Let us consider the scalar soft 
breaking sector of the MSSM: 

- ^soh = mQ__QlQi + mlc u<'*u''i + mic e'=*e'=i + m% d^ldPt + m\ L\Li + m\ hIHi + m% h1H2 

j^j^ j^j^ 11 1 Z 

+ (A?.)LLgjQi + (Ar,)RRU=-n-j + {^.%)KKd%d-j + [^tM^~Lie-j + (Ar,-)LRnL>S- 

+ (Af^)LRdL-dS + --- (1) 

where we have used the standard notation for the MSSM fields and have explicitly written down the 
various A parameters. Notice that the different Alr include the contributions from the trilinear 
terms with the corresponding Higgs vev. Therefore the trilinear couplings in the second line 
represent only the real Higgs couplings and not their vacuum expectation values. 

Consider SU{b) to be the relevant symmetry at the scale where the soft terms are generated. 
Then, taking into account that matter is organized into the SU(5) representations 10 = (g, n^, e'^) 
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Relations at the weak scale 


Relations at AfcuT 






« {ml/ml) (5^^.)rR 




('5.'^-)rR « (^^LL 






/le _ Ad 
^^0 - "^J»o 



TABLE I: Links between various transitions between up-type, down-type quarks and charged leptons for 
SU(5). The suffix '0' implies GUT scale parameters. 

and 5 = (/, d*^), one obtains the following relations 

2 2 2 2 /o\ 

'^Q = ^i- = ^u- = "^10 (2) 

2 2 2 fo\ 

=mL = m^ (3) 

^ij = ^ji ■ (4) 

Eqs. ©-([I]) are matrices in flavor space. These equations lead to relations between the slepton 
and squark flavor violating off-diagonal entries Ajj. These are: 

(Ar,.)LL = (A« )rr = (Af^.)LL = (A',,.)rr (5) 
(Af^.)RR = (Ai,)LL (6) 
(Af,.)LR = (A^,)lr = (A^,.)rl (7) 

These GUT correlations among hadronic and leptonic scalar soft terms are summarized in the 
second column of Table HI Assuming that no new sources of flavor structure are present from 
the SU{5) scale down to the electroweak scale, apart from the usual SM CKM one, one infers 
the relations in the first column of Table [I] at the low-energy scale. A comment is in order when 
looking at Table [H Mass insertions for down-type quarks and leptons in Table [I] always exhibit 
opposite "chiralities" , i.e. LL insertions are related to RR ones and vice- versa. This stems from 
the arrangement of the different fermion chiralities in SU (5) multiplets (as it clearly appears from 
the last column in Table H]). This restriction can easily be overcome if we move from SU{5) to 
left-right symmetric unified models like SO(IO) or the Pati-Salam (PS) case, with gauge group 
SU{4), X SU{2)l X SU{2)r. 
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III. WEAK-SCALE QUARK-LEPTON RELATIONS 

The exact equality of the mass and trihnear matrices of the different components of a GUT 
multiplet is only true when the GUT symmetry is valid. At low scales, where we perform our 
experiments, these relations are modified. In this section, we analyze how these relations are 
modified after the breaking of the GUT symmetry and the RGE running from Mgut to Myy. 

A. Definition of the SCKM basis 

In low-energy phenomenology we obtain the different mass insertion bounds in the so-called 
super-CKM (SCKM) basis, i.e. the basis where through a rotation of the whole superfield (fermion 
-|- sfermion) we obtain diagonal Yukawa couplings for the corresponding fermion field [13]. At low 
energies the SCKM bases for quarks and leptons are unrelated. However in a GUT theory quarks 
and leptons unify and thus we would expect the quark and lepton Yukawa couplings to be equal 
at the GUT scale and hence their SCKM bases to unify. Unfortunately this is not true when we 
evolve the Yukawa couplings using the SM (or MSSM) Renormalization Group Equations (RGEs). 
This discrepancy, that could be solved with some GUT breaking effects or new non-renormalizable 
contributions, implies a relative rotation between the quark and leptonic SCKM basis. While in 
principle this is an obstacle for a model-independent analysis, we argue below that its effects are 
not expected to change the order of magnitude of our bounds. 

Let us consider for example the case in which we obtain a MI bound from a leptonic process 
at low energies. This bound is given in the basis of diagonal charged lepton masses. The next 
step in our scheme consists in using the RGE equations to evolve this bound to the GUT scale. 
The resultant bound is a bound on the element in the sfermion mass matrix common for the 
GUT multiplet including both quark and leptons at Mqut in the basis of diagonal charged lepton 
Yukawa matrices} However, if we want to convert this bound into a bound on the corresponding 
squark MI we must take into account a possible misalignment between the charged-lepton and the 
down-quark Yukawa matrices just below the GUT scale. GUT symmetry breaking effects could 
introduce different corrections to the quark and lepton Yukawa matrices. For example, it is well- 
known that the minimal SU (5) Lagrangian with a 5 and a 5 Higgs representations is not realistic 
as the first two generations of the down quark and charged lepton fermion masses do not follow 

^ An additional (small) RGE-induced rotation can be necessary in the presence of neutrino Yukawa couplings. We 
neglect its effect in our discussion. 
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the GUT relations. This means that quark and lepton Yukawa matrices are not the same below 
Mqut and they are not simultaneously diagonalizable. Taking into account these facts, we must 
rotate our bounds from the leptonic sector to the basis of diagonal quark Yukawa matrices before 
we compare them with the bounds obtained from quark processes at low energies. 

In a complete GUT model we would be able to obtain both the quark and lepton Yukawa 
matrices after the breaking of the GUT symmetry. In such a theory it would be possible to 
compute the relative rotation between the basis of diagonal quark and the basis of diagonal lepton 
Yukawa matrices exactly. On more general grounds, as we do not have a complete GUT theory yet, 
we have a new uncertainty introduced in the translation of the bounds between quark and leptonic 
Mis, due to this relative rotation between the leptonic and quark Yukawa couplings. To see the 
effect of this 'new' rotation matrix, which can be represented by a unitary matrix, U(*'), while 
translating the bounds from one sector to another sector (leptonic to hadronic or vice- versa), we 
consider one particular case. Let ^l^ij, of any 'chirality', have a bound at the GUT scale, given by 
^rnax • This clement A -.^^ corresponds to a combination of matrix elements in the basis of diagonal 
quark Yukawas, , given as 

i^?.i = iE^if l<A™^ (8) 

k,l 

where U^"?') would represent the rotation between the basis of diagonal quark and lepton Yukawas. 
Given that A*^-^) and U('J^) have different origins, we do not expect large cancellations between 
different terms in this sum and therefore, "barring accidental cancellations", we can apply the 
bound on AfJ^ to each of the individual elements in this sum 

< A"""^, (9) 

Once more, in principle it is necessary to know the relative rotation u(^') to apply the bound on 
A[^f . In SU(5) GUT, requiring that the successful b-T unification is not accidental, we can expect 
that \V^f\ ~ 1 and \V^^-\Av}^^\ < 1. The relat ive misalignment for the first two generations 
is much more uncertain. However, we consider reasonable to expect that the different families 
are approximately arranged in the usual way and (c, s) quarks are mostly associated with (fj,, u^j) 
and {u,d) with {e,Ue). So, we make the "reasonable assumption" that diagonal elements in these 
rotations, V^^''\ are 0{1), where 0{1) may mean l/\/2 or even 1/2 but not much smaller. Therefore, 
under these conditions we have that 

< A'""^, (10) 



(¥)* aw) t/w) 
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and thus, order of magnitude MI bounds in the leptonic sector can be safely translated into order 
of magnitude MI bounds in the quark sector, taking into account that these bounds are not precise 
within factors of two. 

In the following we present a realistic example to show the possible effects of these rotations. 



A popular solution to the first two 
by the Georgi-Jarlskog mechanism 



generations Yukawa unification problem in SU(5) is provided 



14l ]. There one introduces an additional Higgs representation 



which is 45 dimensional contributing mainly to the second generation Yukawa couplings. This 
Higgs gets a vev along the {B — L) direction and therefore breaks the charged lepton-down quark 
symmetry. Assuming simple symmetric Yukawa textures, the down quark and charged lepton 
Yukawa matrices are given as follows: 



/ A3 A3\ 



Mgut 



oc 



A3 A2 A2 
A3 A2 1 



OC 



A3 A3 \ 
A3 3A2 A2 

A3 A2 1 y 



(11) 



where A is a parameter of the order of the Cabibbo angle. In this way, we preserve the successful 6— 
r unification and correct the bad GUT relations for the second and first generations. Clearly, if the 
Yukawa matrices have this structure, the basis of diagonal down-quark Yukawas is different from the 
basis of diagonal charged-lepton Yukawas. Therefore the rotation due to the misalignment 

between quarks and leptons is non-trivial. However, it is very easy to check that up to terms of 
0(A3) this rotation is given by 



I 



V 



1 



2A2 
2A 





2A 



2A" n 

1 



(12) 



Using this rotation and barring accidental cancellations in Eq. ([9]), it is evident that Eq. (jlOp is an 
excellent approximation in this case and we can safely relate off-diagonal elements in the squark 
and slepton mass matrices. 



B. RG evolution 



The SU(5) GUT relations between squark and slepton masses and trilinear couplings in Eqs. dJ])- 
([7|) are only valid at the SU(5) scale. As it is well-known, one could expect the long RG running 
from the GUT scale down to the weak scale to modify these relations significantly. However, the 
crucial aspect is that the flavor violating entries (i 7^ j) are not significantly modified due to this 
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running. In fact, even in the presence of right handed neutrinos in a seesaw mechanism, this 
statement remains true up to factors order one. This approximate non-modification of the Aj^j's, 
allows us to recast them at the weak scale in terms of Si^j. In the present section, we will try 
to quantify our statements using semi-numerical solutions of the RGE within the SU(5) model we 
have been discussing so far. We will consider the effects of adding right-handed neutrinos a bit 
later on. 

In the RG evolution of diagonal elements of sfermion mass matrices we neglect small Yukawa 
couplings and keep only Y^, and Yt. Then, for the first two generations of squarks and sleptons. 



the following simple expressions hold in terms of the high scale soft parameters, m^Q 



mi 



ml 



and i^^^ ^ general expression with different m^Q and m| see Table I in 

{ml)i^2{Mw) ^ ml + 6.5 M^/^ {-mlhAMw) ^ + 6.1 Mf^^ 
{ml)i^2{Mw) ^ ml + 0.15 M^^^ {ml)i,2{Mw) ^ m^ + 0.5 M^^^ 
(m^)i,2(Mv 



ml + 6.2 Ml 1 2 



(13) 



Third generation masses receive contributions proportional to Yukawa couplings which depend on 
the value of tan (3. For low tan /3 ~ 5 we have 

{ml^)^{Mw) ^ 0.6 ml + 5.5 mI/2 {ml)3{Mw) ^ ml + 5.8 mI/2 



{ml)3{Mw) ^ ml + 0.16 MI/2 (mlUMw) ^ + 0.5 MI/2 



while for tan /3 ~ 30 we have 

(mlUMw) 
(mlUMw) 



0.2 ml + 4.1 Ml/2 ' 



(14) 



0.5 ml + 5.2 Ml/2 
0.8 ml + 0.12 mI/2 



(mlUMw) ^ 0.8 ml + 5.14 mI/2 
{ml)3{Mw) ^ 0.92ml + 0.5MI/2 



{ml)s{Mw) ^ 0.2 + 4.2 Ml/2 



(15) 



RGE evolution for off-diagonal elements is solely dependent on Yukawa couplings and does not 
include gaugino contributions. Now, if we neglect first two generations Yukawa couplings, in the 

basis of diagonal down-quark Yukawa couplings we have: 

d{mjj^E)i^j 



(167r2 
(IGvr' 



dt 








' dt 




.2.d{ml] 





-2 {ml^E)3j - 2 {ml^E)',3 6,3 
- Yb {ml)-i, 6,3 + Y,^ (m|),3 6,3 + Y,^ V,*Vtk (m^kj + 
yt'VtlVtjimDik +2Y;'V^Vt, [{ml)33 + mj,^) + 2 [y^Y^^ 

-Y^{ml)3j6,3 - Y^{ml),35,3 



(16) 
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where t = logMcuT/Q- We are interested in two different aspects of these equations. First we 
want to determine how a non- vanishing off-diagonal entry in the sfermion mass matrices at Mqut 
is modified due to RGE evolution from Mqut to M\y. Then, we are also interested in the size 
of the off-diagonal entries generated by the running from Mqut to Mw in the case of exactly 
vanishing off-diagonal entries at Mqut- 

Regarding the evolution of GUT elements, we can neglect the non-diagonal CKM elements and 
then, we see from Eq. (jl6p that the (A-^)i2,2i are not modified at the leading order. Furthermore, 
in the low tan/3 limit, off-diagonal elements in m^, or m'^ in the SCKM basis are not significantly 
changed through RG evolution. Thus, one can safely assume Eq. ([6]) to be valid at any scale at 
low tan/3. For large tan /3 ~ 40, the effect of Y^^ is only relatively important and in the leading log 
approximation we obtain 

{ml),=s^,{Mw) ^ (mi),=3^,(Mp^) Exp. (^--L^(2Y,' - Y^) log ^ 0.84 {mi)i=^^,{Mw) 

(17) 

and therefore it can account at most for a (15-20)% relative change. Something analogous happens 
in the left-handed squark mass matrix. For low tan /3 the effect of the top Yukawa coupling in the 
third row elements would be 

{ml),=^^^{Mw) ^ (m|)i=3^,(MH.)Exp. (^--i-^y/ log ^^l) ^ 0.86(m|)i=3^,(MH/), (18) 

In the large tan (3 region, we have to multiply the argument in the exponential by a factor {Y^ + 
Y^ — 2Y^)/Yf^, which for tan/3 = 40 amounts only to a factor of 1.16. Therefore, the relative 
change between (mQ)j=3^j and (m|;)j=3^j is practically the same as the previous case. 

A second feature of these RGE equations is that, even if we start from a flavor blind situation 
at the GUT (or Plank) scale, flavor off-diagonal entries in the LL squark mass matrix at the weak 
scale (via the CKM matrix) are unavoidable. The running from Mqut to the weak scale M\y gives 
rise to the following FV effects 

If we take Mqut — 2 x 10-*^^ GeV, It — 1, — ^^^1/2 + "^0 the hmit M1/2 ^ mo ~ aoi we 

have ^ —0.2 V^jYtj. More generally we can say that {^'11)^ — — c V^*Vtj with c a numerical 

coefficient depending on (mo,Mi/2,ao) taking values between 0.1 and 1. As we will discuss in 
following sections, these RGE induced FV effects have an important phenomenological impact, 
specially when FV entries in the RR squark sector are also present. This is due to the fact that. 
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as we will see, /S.Mb,k mass differences are much more sensitive to contributions proportional to 

(^LL X (^RR- 

Finally, for the last equation in ([7]) , the RGE equations for off-diagonal down quark and charged 
lepton trilinear couplings are 

(levr^)^^^!^ = (f <7i + 35i + loi) {Y,% - {Y,%, - 4^^ (y/)^3 s^, (20) 



dt 



(16^2):n^e^ = (3gl + Sgj) {Yf},, - 5Y,HK%j - 4^,^ (^/k ^,3 - if (^^^ 



-3Y^ {Yf-)ij. 



e 



From the equations above, it's clear that the trilinear couplings scale almost in the same manner 
as the corresponding Yukawa couplings. Thus the corresponding relation given in Table I holds 
well especially in the limit of small tan/3 ~ 25, where we can neglect the bottom and r Yukawas. 
For large tan/3 and for 13 and 23 Mi's, we expect similar corrections ~ (15 — 20)% as we have 
mentioned above. Given that the MI bounds are intrinsically order of magnitude bounds, we can 
safely consider these off-diagonal elements to evolve analogously to the Yukawa couplings from M]y 
to Mgut- 

Let us now briefly mention the implication of adding right-handed neutrinos as in a seesaw 



mechanism, 
literature Q, 



?he RG effects and subsequent flavour violation have been studied in detail in the 
17l |. The implication of these effects on the relations Eqs. ([5][7I) depends crucially on 



the strength of the neutrino Dirac yukawa couplings Y^ and the "mixing" they carry. For example, 
a rough estimate of the weak scale (A^^)jj is given by 

{m^f,^^{Mw) ^ {m^%^{MGUT) - ^ J^i^.Y^hiSml + A^) log (21) 

where Yi, represents the neutrino Dirac Yukawa couplings. From the equation above we see that 
RG effects could "generate" large As if the mixing in Y^ is large. In the general case, it is very 
difficult to separate the RGE-induced As from the RGE-modified GUT off-diagonal elements in 
m| and Yj^. Changes in the original GUT off-diagonal entries and the RGE-induced As depend on 
the neutrino Yukawa couplings. These neutrino Yukawa matrices are present in the RGE between 
the scales Mgut and the mass of the corresponding right handed neutrino. Typically one of the 
eigenvalues in the neutrino Yukawa matrices can be of the same order as the top Yukawa coupling,^ 



but this large eigenvalue is necessarily associated with the heaviest right-handed neutrino 



18|, 



We do not consider here the possibility of having several large eigenvalues in 
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and therefore it decouples early. The main problem regarding the neutrino Yukawa couplings is 
that the mixing diagonalizing Y,yYj in the basis of diagonal Yf,Y^ is unknown and could be large. 
Nevertheless, we can discuss some cases. 

The simplest situation corresponds to the case where the mixing diagonalizing the charged 
lepton and neutrino Yukawa matrices are both small. In this case, the large neutrino mixings are 
generated through the seesaw mechanism and they play no role in the RGE evolution of slepton 
matrices. The RGE evolution of m? is then very similar to the previous cases without the right 
handed neutrinos. In the low tan f3 region we can completely neglect the effects of Y(,Y^ in the 
RGEs. Then in the basis of diagonal neutrino Yukawas (approximately equal to the basis of 
diagonal charged-lepton Yukawas) the off-diagonal elements would be 

{ml),=sMMw) ^ iml),=sMMw) Exp. (^-i^(y^2 ^ ^ iml),=3^j{Mw), (22) 

taking Mr^ as low as 10^° GeV. 

In the large tan f3 region, the effects of Eq. ()22p and Eq. (I17p have different signs and partially 
cancel out. So, we can also expect these off-diagonal elements to remain of similar magnitude 
between Mgut and Mw. Notice also that the (1,2) mass insertions remain unchanged if we have 
only one large neutrino- Yukawa coupling corresponding to the third generation. 

The second situation corresponds to large mixings in the neutrino Yukawa matrices in the basis 
of diagonal charged-lepton Yukawas. This case depends strongly on the mixings and size of the 
neutrino Yukawa couplings. There are two simultaneous effects, the creation of a new MI due to the 
RGE and the change from RG evolution of the GUT delta. These effects cannot be sorted out in 
general and thus we cannot describe the evolution of the GUT delta without specifying completely 
the Yukawa matrices. In this situation, it is difficult to correlate quark and leptonic deltas. Finally, 
the effects of neutrino Yukawas on the relations between charged-lepton A-parameters is expected 
to be similarly small even for large enough (top-quark like) Yukawa couplings. 

We conclude noting that, in general, the RG effects do not significantly modify the GUT rela- 
tions for off-diagonal elements already present at Mqut- However in the presence of right-handed 
neutrinos and large mixings in the neutrino Yukawa couplings, the new contributions to the off- 
diagonal entries in the slepton mass matrices can destroy the quark-lepton correlations and the 
analysis becomes model dependent. 
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IV. CONSTRAINTS 

In this section, we list the constraints we have imposed on the SUSY parameter space before 
starting the analysis of flavor physics in the leptonic and hadronic sectors. 

A. Direct SUSY search 

The model we use in this work is a minimal departure from the usual Constrained MSSM 
(CMSSM) where we add a single flavor changing entry in the sfermion couplings at Mqut- The 
usual CMSSM contains (assuming vanishing flavor blind SUSY phases) five parameters: 
mo, tan/3 and sign(^)).'^ In addition we include the single MI relevant to the process(-es) of 
interest at the GUT scale. Here, we take the /i sign positive as required by the 6 — > 57 and muon 
anomalous magnetic moment constraints. Therefore the SUSY parameters that enter our analysis 
for a given MI are only (mo,mi/2i^Oitan/3,A^). We scan the values of these parameters in the 
following r.n,.s: M,., < m GeV, ™„ < 380 GeV, | A„| < 3™„ and 5 < fan ^ < 15^ The bound on 
the ^0 parameters is set to avoid charge and / or colour breaking minima [20|] ■ This would typically 
correspond to the following highest mass scales: slepton masses as large as ~ 400GeV and 
squark masses as large as niq ~ 550 GeV; both possibly observable at LHC. 

At the low scale, we impose the following constraints on each point: 



Lower bound on the light and pseudo-scalar Higgs masses 



2l|; 



The LEP constraints on the lightest chargino and sfermion masses 



The LEP and Tevatron constraints on squark and gluino masses 22 1. 



We consider the MSSM with conserved R-parity and thus the lightest neutralino (that coincides 
with the lightest supersymmetric particle - LSP) provides an excellent dark matter candidate. 

B. & — !■ S7 

B{B ^sT) is particularly sensitive to possible non-standard contributions and it provides a 
non-trivial constraint on the SUSY mass spectrum given its precise experimental determination 
and the very accurate SM calculation at the NNLO 23|]. When witnessing such a light SUSY 



^ Notice, we are not enforcing here the minimal supergravity relation between the A and B parameters, B — A — 1 
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spectrum (we remind that we take Mi 12 < 160 GeV, mo < 380 GeV) a legitimate worry is 
whether the bounds on B{B Xg^) are respected. In this work we choose /i > which impHes 
destructive interference between chargino and charged Higgs contributions (and is also preferred by 
the {g — 2)^ constraints). We have explicitly checked that the combined chargino and charged Higgs 
contributions satisfy the B{B Xg^) constraints. For these points, we checked simultaneously 
that gluino contributions satisfy by themselves the B{B Xg'j) constraints. Then, these gluino 
contributions set a bound on the 6'^^ Mis. Notice that, in this way, we are not allowing the possibility 
of an accidental cancellation of charged-Higgs and chargino contributions with gluino ones. 



C. The lightest Higgs boson mass 

The non-observation of the lightest neutral Higgs boson (h^) at present colliders is already a 
stringent constraint on the MSSM parameter space 211]. Even if the mass depends on the whole 
set of MSSM parameters (after the inclusions of loop corrections), m^o mainly depends on (and, 
indeed, grows with) the left-right mixing term in the stop mass matrix M/"-^ = mt{Au — n/ tan (3), 
on the average stop mass Mg and on tan/3. In particular, it is well known that values of tan/3 < 2 
are strongly disfavored. Taking into account that the lower bound on the /i" mass changes with 
the SUSY parameters, we have explicitly checked that the experimental limits 
our parameter space. 



21I ] are fulfilled in 



D. (5-2)^ 



The possibility that the anomalous magnetic moment of the muon [a^ = (g — 2)^/2], which has 



been measured very precisely in the last few years 



24| , provides a first hint of physics beyond the 



SM has been widely discussed in the recent literature. Despite substantial progress both on the 



experimental and on the theoretical sides, the situation is not completely clear yet (see Ref. 25 1 
for an updated discussion). 



Most recent analyses converge towards a 2a discrepancy in the 10 ^ range 



exp 



(2± 1) X 10" 



(23) 



The main SUSY contribution to usually provided by the loop exchange of charginos and 

sneutrinos. The basic features of the super symmetric contribution to are correctly reproduced 
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by the following approximate expression: 



-1.5fi;^ii?)f!^V('i^if!^ , (24) 



1 X 10-9 ■ V 10 / V m;^ ^ 
which provides a good approximation to the full one- loop result [26 1. 



The most relevant feature of Eqs. (j24p is that the sign of a^^^^ is fixed by the sign of the n 
term so that the fi > region is strongly favored. This is specially true for the Standard Model 
prediction which uses the data from e^e^ collisions to compute the hadronic vacuum polarization 
(HVP). This predicts a smaller value than the experimental result by about 3 a. In case one uses 
the r data to compute the HVP, the discrepancy with SM is reduced to about 1 a, but it still 
favors a positive correction and disfavors strongly a sizable negative contribution. Thus, taking 
fj, > 0, the region of parameter space considered in this analysis satisfies the constraint of Eq. (j23p . 

E. Electroweak Precision Observables (EWPO) 

The good agreement between the SM predictions and the electroweak precision observables 
(EWPO) points to a decoupling of new physics contributions to these precision observables. As 
we also consider light superpartners, we need to take into account the tight constraints on the 
supersymmetric spectrum emerging from this agreement. 

Several recent and thorough analyses for the MSSM are available in the literature 0. The 
most relevant effect is due to the mass splitting of the superpartners, and in particular of the third 
generation squarks. Indeed, large splitting between rrif,^ and rht^ would induce large contributions 
to the electroweak p parameter. This universal contribution enters the boson couplings and the 
relation between Mw, and a and is therefore significantly constrained by present data. In the 
cases of the M^-boson mass and of the effective weak mixing angle sin'^ ©W' ^"-"^ example, a doublet 
of heavy squarks would induce shifts proportional to its contribution to p: 6Mw /Mw ~ 0.72 A/9 
and 5sin2 9^^ ^ -0.33 Ap. 

In our analysis, we have required that A^p^'^^ < 1.5 x 10^'^ and we have checked that no relevant 
constraints arise from EWPO, as it is confirmed by the thorough analysis (relative to the CMMSM 



framework) in Ref. 27 1. 



V. MASS INSERTION BOUNDS FROM HADRONIC PROCESSES 



The comparison of several hadronic flavor-changin 
be used to bound the Mis in the different sectors 



processes to their experimental values can 



lot]. In these analyses it is customary to 
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Observable 


Measurement /Bound 


Ref. 






Sector 1- 


-2 






AMk 




(0.0 - 5.3) X 10-3 GeV 


[30] 




e 




(2.232 ±0.007) X 10^^ 


[30] 




\{£'/£)susy\ 




< 2 X 10-2 


- 






Sector 1- 


-3 










(0.507 ± 0.005) ps-i 






sin 2/3 




0.675 ±0.026 






cos 2/3 




> -0.4 


[32] 






Sector 2- 


-3 






{s + d)-/){E^ > 2.0 GeV) 


(3.06 ±0.49) X 10-"^ 






{s + d)j){E^ > 1 


8 GeV) 


(3.51 ±0.43) X lO-'^ 


m 


BR(6 


-> s-f){E^ > 1.9 Gel/) 


(3.34 ±0.18 ±0.48) x lO"'^ 


[35] 








0.004 ±0.036 


[31] 


BR(6 ^ sl+ 


r)(0.04 GeV < q 


2 < 1 GeV) 


(11.34 ±5.96) X lO-'^ 


[36,37] 


BR{b ^ sl+l-){l GeV < g2 


< 6 GeV) 


(15.9 ±4.9) X lO-'^ 


[36,37] 


BR(& ^ sl+l 


-)(14.4 GeV < g2 


< 25 GeV) 


(4.34 ± 1.15) X lO-'^ 


[36,37] 




Acpib ^ sl+l-) 




-0.22 ±0.26 


[30] 




AMb, 




(17.77 ±0.12) ps-i 


[38] 



TABLE II: Measurements and bounds used to constrain the hadronic S s. 

consider only the dominant contributions due to gluino exchange which give a good approximation 
of the full amplitude, barring accidental cancellations. In the same spirit, the bounds are usually 
obtained taking only one non-vanishing MI at a time, neglecting the interference among Mis. 
This procedure is justified a posteriori by observing that the MI bounds have typically a strong 
hierarchy, making the destructive interference among different Mis very unlikely. 

The effective Hamiltonians for AF = 1 and AF = 2 transitions including gluino contributions 
computed in the MI approximation can be found in the literature together with the formulae 
of several observables \^. Even the full NLO calculation is available for the AF = 2 effective 



Hamiltonian 



28|, 



|29|]. 

In our study we use the phenomenological constraints collected in Table [Til We use the same 
set of SUSY parameters described in the previous Section, so that hadronic and leptonic Mis are 
related as discussed in Section [Til In particular: 
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Sector 1-2 The measurements of ^Mk, s and e' /e are used to constrain the ('^12)^^ with 
{A^B) = {L,R). The first two measurements, AMk and e respectively bound the real 
and imaginary part of the product ((^12) (<^i2)- case of AMk, given the uncertainty 

coming from the long-distance contribution, we use the conservative range in Table [Til The 
measurement of e'/e, on the other hand, puts a bound on Im.{6i2)- This bound, however, is 
effective in the case of the LR MI only. Notice that, given the large hadronic uncertainties in 
the SM calculation of e' /e, we use the very loose bound on the SUSY contribution shown in 
Table [III The bounds coming from the combined constraints are shown in Table IIIII Notice 
that, here and in the other sectors, the bound on the RR MI is obtained in the presence of 
the radiatively-induced LL MI given in Eq. (|19p . The product [S12) (^12) generates 
left-right operators that are enhanced both by the QCD evolution and by the matrix element 
(for kaons only). Therefore, the bounds on RR Mis are more stringent than the ones on LL 
Mis. 

Sector 1-3 The measurements of AMg^ and 2(5 respectively constrain the modulus and the 
phase of the mixing amplitude bounding the products ((^is) ("^is) ■ -^"-"^ sake of simplicity, 
in Table Hill we show the bounds on the modulus of (^^is) only. 

Sector 2-3 This sector enjoys the largest number of constraints. The recent measurement 
of AMb^ constrains the modulus of the mixing amplitude, thus bounding the products 
I ('^23) (^23) I- Additional strong constraints come from /S.B = 1 branching ratios, such as 
b —>■ s'y and b — > sl^l~ . Also for this sector, we present the bounds on the modulus of (^(^fs) 
in Table HUl 

All the bounds in Table Hill have been obtained using the NLO expressions for SM contributions 



and for SUSY where availab 
lattice calculations 



391, 



401, 



41 



e. Hadronic matrix elements of AF = 2 operators are taken from 
□ 

42l | . The values of the CKM parameters p andn are taken from the 



431]. This conservative 



UT/ii analysis in the presence of arbitrary loop-mediated NP contributions 
choice allows us to decouple the determination of SUSY parameters from the CKM matrix. For 
6 — > 57 we use NLO expressions with the value of the charm quark mass suggested by the recent 
NNLO calculation 



23]. For the chromomagnetic contribution to e' /e we have used the matrix 
element as estimated in Ref. 44]. The 95% probability bounds are computed using the statistical 
method described in Refs. 3, 451 • 

Concerning the dependence on the SUSY parameters, the bounds mainly depend on the gluino 
mass and on the "average squark mass" . A mild dependence on tan (3 is introduced by the presence 
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tj\AB 


LL 


LR 


RL 


RR 


12 


1.4 X 10-2 


9.0 X 10-^ 


9.0 X 10-5 


9.0 X 10-3 


13 


9.0 X 10-2 


1.7 X 10-2 


1.7 X 10-2 


7.0 X 10-2 


23 


1.6 X 10-1 


4.5 X 10-3 


6.0 X 10-3 


2.2 X 10-1 



TABLE III: 95% probability bounds on | {^ij) I obtained using the data set described in Section ITVl See 
the text for details. 



Process 


Present Bounds 


Expected Future Bounds 


BR(/i^ 6 7) 


1.2 X 10-11 


0(10-13-10-14) 


BR(/^ eee) 


1.1 X 10-12 


0(10-13-10-14) 


BR(^ ^ e in Nuclei (Ti)) 


1.1 X 10-12 


0(10-18) 


BR(t -^6 7) 


1.1 X 10-^ 


0(10-8) 


BR(r ^ eee) 


2.7 X lO-'^ 


0(10-8) 


BR(t —ten /i) 


2. X lO-'^ 


0(10-8) 


BR(t ^ /i7) 


6.8 X 10-^ 


0(10-8) 


BR(t M M 


2 X lO-'^ 


0(10-8) 


BR(t —* fiee) 


2.4 X 10-^ 


0(10-8) 



TABLE IV: Present and Upcoming experimental limits on various leptonic processes at 90% C.L. 

of double Mis {^fj) ('^■^■')l_r chromomagnetic operators. This dependence however becomes 
sizable only for very large values of tan (3. 



VI. MASS INSERTION BOUNDS FROM LEPTONIC PROCESSES 

In this section, we study the constraints on slepton mass matrices in low energy SUSY imposed 
by several LFV transitions, namely li — > Ij'j, li Ijhh and fi-e transitions in nuclei 



46(]. The 



present and projected bounds on these processes are summarized in Table IIVI These processes 
are mediated by chargino and neutralino loops and therefore they depend on all the parameters 
entering chargino and neutralino mass matrices. In order to constrain the leptonic Mis 6^^, we 
will first obtain the spectrum at the weak scale for our SU(5) GUT theory as has been mentioned 
in detail in section IIVI Furthermore, we take all the flavor off-diagonal entries in the slepton 
mass matrices equal to zero except for the entry corresponding to the MI we want to bound. To 
calculate the branching ratios of the different processes, we work in the mass eigenstates basis 
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Type of 




fi ^ eee 


M - 


-> e conversion in Ti 


LL 


6 X 10-^ 


2 X 10"^ 




2 X 10-3 


RR 




0.09 






LR/RL 


1 X lO"'"^ 


3.5 X 10-5 




3.5 X 10-5 



TABLE V: Bounds on leptonic from various /i — > 6 processes. The bounds are obtained by making a scan 
of mo and Afi/2 over the ranges mo < 380 GeV and Afi/2 <160 GeV and varying tan (3 within 5 < tan /? < 15. 
The bounds are rather insensitive to the sign of the /i mass term. 

through a full diagonalization of the slepton mass matrix. So, imposing that the contribution of 
each flavor off-diagonal entry to the rates of the above processes does not exceed (in absolute value) 
the experimental bounds, we obtain the limits on the (5*-^'s, barring accidental cancellations. 

The process that sets the most stringent bounds is the li Ijj decay, whose amplitude has the 
form 

T = mi^e^ujip - q)[iq''axMLPL + ArPr)]u,{p) , (25) 

where p and q are momenta of the leptons 1^ and of the photon respectively, Pr^l = ^(1 =t 75) and 
Al^r are the two possible amplitudes entering the process. The lepton mass factor mi. is associated 
to the chirality flip present in this transition. In a super symmetric framework, we can implement 
the chirality flip in three ways: in the external fermion line (as in the SM with massive neutrinos), 
at the vertex through a higgsino Yukawa coupling or in the internal gaugino line together with a 
chirality change in the sfermion line. The branching ratio of li Ij^ can be written as 



BRjlj Ijj) _ 487r3a ,j ^ , 1 .ij 



2^ 



with the SUSY contribution to each amplitude given by the sum of two terms Al^r = A^ r + _r- 
Here A^ ^ and denote the contributions from the neutralino and chargino loops respectively. 

Even though all our numerical results presented in Tables IVllVIII are obtained performing an 
exact diagonalization of sfermion and gaugino mass matrices, it is more convenient for the discussion 
to use the expressions for the li Ijj amplitudes in the MI approximation. In particular, we treat 
both the slepton mass matrix and the chargino and neutralino mass matrix off-diagonal elements 
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Type of ,513 


r — !■ e 7 


r — > e e e 




LL 


0.15 






RR 








LR/RL 


0.04 


0.5 





TABLE VI: Bounds on leptonic from various r ^ e processes obtained using the same values of SUSY 
parameters as in Table fVl 



as mass insertions.^ In this approximation, we have the following expressions 



LL 



+ 



/in (02 ) + /ic (02 ) + ^ffa f -^Sn (a2 , &) + /2c (a2 , ^) 

fsnjai) f2niai,b) 
m? (/x2-M2) 



+ tan ^14/ /in(ai) + /i-/Vfi tan /3 



47r m? 



f^) 2/2„(ai), 



(26) 



2 



+ 



LR 



4/in(ai) +AiMitan/? 
Ml' 



/3n(ai) _ 2/2n(ai,fc) 



2 hn{ai) 



(27) 



where 0vK is the weak mixing angle, ai^2 = -^^1,2/"^ ; ^ = ^ /mj and fi[c,n){xiy) = /i(c,n) 
fi{c,n)iy)- The loop functions /j are given as 

-17x^ + Ox^ + 9x - 1 + 6x2(x + 3) lnx)/(24(l - xf), 
-5x^ + 4x + 1 + 2x(x + 2) lnx)/(4(l - xf), 
(1 + 9x - + 6x(x + 1) lnx)/(3(l - xf), 

x^ - 9x2 + 9x + 1 + g^^^ _^ x)/(6(l - x)^), 



X — 



/in (2; 
/2n(2; 
/3n(a; 
flc{x 
/2c (2; 



4x + 5 + 2(2x + 1) In x)/(2(l - x)^) . 



(28) 



We note that all {^Ijjj^^ contributions with internal chirality flip are tan /3-enhanced. On the other 
hand, the only term proportional to arises from pure B exchange and it is completely 

independent of tan/3, as can be seen from Eqs. (I26p and (I27p . Therefore the phenomenological 



This approximation is well justified and reproduces the results of the full computation very accurately in a large 
region of the parameter space [4^ . 
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Type of (5^3 


r ^ /Z7 


T ^ fl fl fl 


T ^ flee 


LL 


0.12 






RR 








LR/RL 


0.03 




0.5 



TABLE VIL Bounds on leptonic from various t ^ fi processes obtained using the same values of SUSY 
parameters as in Table fVl 

bounds on (^L) depend on tan/3 to some extent, while those on (^L) do not. The bounds 

V / LL \ / LPt 

on LL and RL Mis are expected to approximately fulfill the relation 



(sl ) ,,!^tan/3 (4.) 



This is confirmed by our numerical study. 

The sector requires some care because of the presence of cancellations among different 
contributions to the amplitudes in regions of the parameter space. The origin of these cancellations 
is the destructive interference between the dominant contributions coming from the B (with internal 
chirality flip and a flavor-conserving LR mass insertion) and BH^ exchange Q, I^]. We can see 



this in the MI approximation if we compare the tan /3 enhanced terms in the second line of Eq. (j26ll 
with the tan/3 enhanced terms in Eq. (|27p . Here the loop function fsiai) corresponds to the pure 
B contribution while f2n{o.i,b) represents the BH^ exchange. These contributions have different 
relative signs in Eq. (f26ll and Eq. (f27|) due to the opposite sign in the hypercharge of SU (2) doublets 
and singlets. Thus, the decay k — > Ijj does not allow to put a bound on the RR sector. We can 
still take into account other LEV processes such as li Ijhh and fi-e in nuclei. These processes 
get contributions not only from penguin diagrams (with both photon and Z-boson exchange) but 
also from box diagrams. Still the contribution of dipole operators, being also tan/3-enhanced, is 
dominant. Disregarding other contributions, one finds the relations 



2 



log- 



t2 

k 



Br{li^lj^) 37r y 
Br{fi - e in Ti) ~ aeBR{pL 67) , (29) 

which clearly shows that li — > /j7 is the strongest constraint and gives the more stringent bounds on 
the different 5ij^s. As we have mentioned above, however, in the case of 6\^^ the dominant dipole 
contributions interfere destructively in regions of parameters, so that Br[li Ij'j) is strongly 
suppressed while Br{fi — e in nuclei) and Br{li Ijlkh) are dominated by monopole penguin 
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(both 7* and Z-mediated) and box diagrams. The formulae for these contributions can be found 



in Ref. 4a|. However, given that non-dipole contributions are typically much smaller than dipole 
ones outside the cancellation region, it follows that the bound on S^^^ from fx eee are expected 
to be less stringent than the one on from /i ^ 67 by a factor ^/oJJStt) l/tan/3 ~ 0.02/ tan /3, 
if the experimental upper bounds on the two BRs were the same. This is partly compensated 
by the fact that the present experimental upper bound on the BR{fj, — > eee) is one order of 
magnitude smaller than that on BR(fj, — > 67), as shown in Tab. IIV[ On the other hand, the 
process BR{fi — e in nuclei) suffers from cancellations through the interference of dipole and non- 
dipole amplitudes as well. These cancellations prevent us from getting a bound in the RR sector 
from the /i-e conversion in nuclei now as well as in the future when their experimental sensitivity 
will be improved. However, the — > 67 and ^u-e in nuclei amplitudes exhibit cancellations in 
different regions of the parameter space so that the combined use of these two constraints produces 



a competitive or even stronger bound than the one we get from BR(fi — > eee) alone [46 1 . 

We summarize the different leptonic bounds in tables IVllVIIi All these bounds are obtained 
making a scan of and Mi^2 o'^^r the ranges mo <380GeV and Mi^2 <160GeV and therefore 
correspond to the heaviest possible sfermions. As expected, the strongest bounds for ^^'^ ^lr 
come always from // — > e7, t ^ fij and t ^ processes. In the case of 5\^^ we can only obtain a 
mild bound for (^i2)j^j^ from /i eee and there are no bounds for (^(523) ^^^^ ^-iid ('^isjj^j^- Notice, 
however, that does not mean that these LFV processes are not effective to constrain the SUSY 
parameter space in the presence of RR Mis. For most of the values of rriQ and Mi 12 there is no 
cancellation and the values of these MI are required to be of the order of the LL bounds. Only 
for those values of mo and Mi /2 satisfying the cancellation conditions a large value of the RR MI 
is allowed. Therefore, we must check individually these constraints for fixed values of the SUSY 
parameters. 



VII. QUARK-LEPTON MI RELATIONS IN GUT SCENARIOS 

In the previous two sections we have collected the MI bounds obtained from the hadronic and 
leptonic processes. In the present section, let us consider a GUT theory, with the corresponding 
GUT symmetric relations holding at the GUT scale. We will focus on the SU(5) case as summarized 
by the relations in Table I. To make a comparison between leptonic and hadronic MI bounds we have 
to take into account that these bounds have a different dependence on the low-energy parameters 
of the theory. On the one hand, the hadronic processes are dominated by gluino contributions and 
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therefore they mainly depend on the gluino mass and the "average squark mass".^ 

On the contrary, as we saw in section IVTl the leptonic bounds depend basically on three param- 
eters: gaugino mass, "average slepton mass" and tan In our model squark and slepton masses 
originate from a common scalar mass mo at the GUT scale and therefore we can relate the average 
squark and slepton masses. 

As we have discussed in section llIIB( the off-diagonal elements of sfermion mass matrices are not 
significantly modified in the RGE evolution from Mqut to M^. However, under some conditions, 
as for example in the presence of large neutrino Yukawa couplings, the RG evolution can generate 
sizable off-diagonal elements in the slepton mass matrices even starting from a vanishing value at 
-^GUT- Clearly these effects are never present in the squark mass matrices, thus breaking the GUT 
symmetric relations. This implies that, given our ignorance on the structure of neutrino Yukawa 
couplings, we have to be careful when applying the MI bounds obtained from quarks to leptons or 
vice- versa. 

In fact, if we obtain a bound on a 6lj MI from a leptonic process at low scales, we can say that, 
barring accidental cancellations, this bound applies both to the mass insertions already present at 
-^GUT and to the mass insertions generated radiatively between Mqut and Mi,^. Therefore we 
can translate this low-scale bound into a bound on the MI at the GUT scale. This bound applies 
also to the squark MI at Mqut and using RGEs we can transport this bound to the electroweak 
scale. For example, in SU{5), we find: 

|(4)rrI < (30) 

The situation is different if we try to translate the bound from quark to lepton Mis. An hadronic 
MI bound at low energy leads, after RGE evolution, to a bound on the corresponding grand- 
unified MI at Mgut, applying both to slepton and squark mass matrices. However, if the neutrino 
Yukawa couplings have sizable off-diagonal entries, the RGE running from Mgut to M]y could 
still generate a new contribution to the slepton MI that exceeds this GUT bound. Therefore 
hadronic bounds cannot be translated to leptons unless we make some additional assumptions on 
the neutrino Yukawa matrices. 

On general grounds, given that SM contributions in the lepton sector are absent and that the 
branching ratios of leptonic processes constrain only the modulus of the Mis, it turns out that all 

the MI bounds arising from the lepton sector are circles in the Re ( 6^, ] -Im ( 6^ ] plane and 
^ ^ V 'J/AB V y/AB 

^ Note that the tan/3 dependence seeps in once we consider the double Mis (^fj)^^ i^^i) lr rl' 
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FIG. 1: Allowed region in the Re((523)Lp{^-Im(<^23)LR pl^^^c using hadronic constraints only (upper left), 
leptonic constraints only (upper right), all constraints (lower left) and all constraints with improved leptonic 
bounds (lower right). 

are centered at the origin. In some cases, the hadronic bounds from B physics constraints are too 
loose and, in principle, this would allow the presence of Mis larger than one. The Mass Insertion 
approximation cannot be trusted when the bounds on the 5 approach values 0(1). Therefore, in 
our analysis we always consider MI values smaller than one.^ In the following we will analyze 

^ Even in the case we consider the possibility of 0(1) Ss, the requirement of absence of tachyonic scalar masses in 
the slepton sector, i.e. (^Si)^^ ^ li provides a bound on the squark iMIs through the GUT-symmetric relations 
among leptonic and hadronic iMIs (see Table |l| 
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the effect of leptonic bounds on the quark mass insertions. For instance, if we had a (A23)^^ 

at the GUT scale, this would have effects both in the r /X7 and b ^ decays. Neglecting 

the effects of neutrino Yukawas that, if present, could generate an additional (^62^^^^ in the RGE 

evolution, and using [623] ~ {mh/mr) (m?/m?) [62^] , a bound on ((523) from the t ^ fj.j 
\ ' LR ^ ^ \ / RL \ / RL 

decay translates into a bound on (A23)^^ thus, into a bound on the SUSY contributions to 



BR(i? ^ Xs7). Similarly, the bound on ((523) would translate into an upper bound for the 

\ / LR 

r — > |U7 branching ratio. 

We present the effect of this GUT correlation in our numerical analysis in Fig. [TJ In the top row, 
we show the allowed region in the Re ((523) -Im((523) plane (larger boxes correspond to higher 

\ / LR \ / LR 

probability densities), using hadronic (left) or leptonic (right) constraints only. We see that the 
present leptonic bounds have no effect on the {^2'i^^^ couplings. This is due both to the existence 
of strong hadronic bounds from 6 — > 57 and CP asymmetries and to the relatively weak leptonic 
bounds here. Even assuming a future bound on BR(r ^7) at the level of 10"*^, attainable at 
B factories, leptonic bounds would marginally improve the hadronic constraints. We remind the 
reader that the LR bounds are basically independent of tan /? and hence this fact does not change 
for different tan/3 values. In Fig. [2] we present the results of the same analysis for ('^23)j^j^- While 
the leptonic bounds do not change with respect to the previous case, the hadronic ones are different 
as SM and SUSY 6 — > 57 and h — > si^i~ amplitudes interfere in the LR case only. 

Similarly, if we have a (A23)^^, the corresponding Mis at the electroweak scale are ('^23)^^^^ 
and ( (523 ) that contribute to AMb^ and r — > fi'y respectively (the impact of ( A23 ) on 6 — > 57 

V / LL ^ / RR 

and b si'^i^ is not relevant because of the absence of interference between SUSY and SM 

contributions) . In Fig. [3] we present the allowed values of Re ( (5^3 ) and Im ( (5^3 ) with the 

V / RR V / RR 

different constraints. The leptonic constraints are quite effective as the bound on the BR{t 
fij) from B-factories is already very stringent, while the recent measurement of AMb^ is less 
constraining. The plots correspond to 5 < tan/3 < 15, thus, the absolute bound on ( 

^23 ) is set 

by tan/3 = 5 and it scales with tan/3 as (^23^^^ ~ (5/ tan /?). 

In Fig. S] we show the results of our analysis for ( 62-3 ) . In this case, there is no appreciable 
improvement from the inclusion of leptonic constraints. In fact, we remind that r — > ^7 is not 
effective to constrain ((523^^^, i.e. the leptonic MI related to (^^23^^^ in our SUSY-GUTs scheme, 
in large portions of the parameter space because of strong cancellations among amplitudes. 



^ Sizable SUSY contributions to AMbs are still possible from the Higgs sector in the large tan /3 regime both within 
[igl . [50I ] and also beyond @] the Minimal Flavor Violating (MFV) framework. However, for our parameter space, 
the above effects are completely negligible. 
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FIG. 2: Same as Fig.[T]for (<S23)r,l- 



The analysis of the constraints on the different ( 6'^ 



13 



Mis is similar to that of the (623^ Mis. 



In this case, the hadronic constraints come mainly from AA/^^ and the different CP asymmetries 
measured at B-factories. The leptonic bounds are due to the decay r — > 67. We present the 
numerical results in Figs. O [6] and [71 Notice that, in spite of comparable experimental resolutions 

on AMbs and AMb^ , the constraints on ( ] are stronger than those on f S23 ) (see Figs. I3l6p . 

V / iRiR V / R,R, 

The reason is that, in the 13 sector, we can make use also of the constraints from 2/3, in addition to 
those relative to AMb^- Moreover, from Fig. [6] we see that the constraints arising from a combined 
analysis of leptonic and hadronic processes are much more effective then the bounds obtained from 
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FIG. 3: Same as Fig.[I]for the ('^23)rr- 

the hadronic and leptonic processes alone. This is due to the fact that the maximal allowed values 
for the hadronic and leptonic deltas in the upper row of Fig. [6] correspond to different values of 
(mo,-/Wi/2) in the two cases. ^ So, the different mo dependence of AM^^ and r ^ 67 provide the 
explaination of their interplay in constraining ( 1 (as it is clearly shown in the lower plot on 
the left of Fig. [6l). On the other hand, the above interesting interplay is not effective in the AM^^ 

* From Eq. (|19p we see that the maximal value of the radiatively induced 5^;^ corresponds to small M1/2 (small rriq) 
and large mo. The largest allowed value for 5^-^ is set by the minimum value of this radiatively induced ^lL' 
large M1/2 and small mo/Mi/2. On the other hand the maximal allowed values from the leptonic delta correspond 
to large M1/2 and large mo. 
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FIG. 4: Same as Fig. [T] for (6^^) 



LL- 



case due to the better bound on the decay r — > /X7 and the absence of analog 2/3 constraints in the 
23 sector. 

Moreover, the leptonic bounds do not have a sizable impact in the RL and LL cases, as clearly 
shown in Figs. [5] and [71 respectively. The LR case is identical to the RL one. 

Finally we analyze the 1-2 sector. In Fig. [8] we can see the allowed values of Re ('^12) ^^^^ ^^'^ 
Im (^Si2^^^- As before, the plot on the left of the upper row corresponds to the allowed values of 
these parameters from hadronic constraints. The dominant hadronic bound comes from ex which 
however is ineffective along the Re ('^12) ^^^^ I™ ('^i2)j^j^ axes. These directions are eventually 
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FIG. 5: Same as Fig.[T]for {Sfs}^^. 

bounded by the milder AMk constraint. The upper right plot represents the values allowed taking 
into account the limits on the branching ratios of the processes /x — > 67, /x — > eee and /i-e conversion 
in nuclei as per the SU(5) relations between ((^12) and (5i2) • As we saw in the previous 
section, the /x — > 67 decay does not provide a bound to this MI due to the presence of cancellations 
between different contributions. We can only obtain a relatively mild bound, ('^i2)j^j^ — 0.09, if we 

take simultaneously into account all the leptonic processes. However, we see in the lower left plot 

~ 2 

that, once rescaled by the factor this bound is more stringent than AMk, so that it cuts the 
tails along the axes. Using the expected bounds for these decays from the proposed experiments. 
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FIG. 6: Same as Fig.[T]for (<5i3)j^j^. 

we obtain the lower right plot in the figure. In these plots the leptonic constraints come from 
the monopole and box contributions to fi ^ eee and ^-e conversion in nuclei and therefore these 
bounds are independent of tan/3. There is a modest improvement of the bounds on (^Si2^^^ which 
however do not take into account possible improvements of the ex constraint. 

In Fig. [9] we present the allowed values of Re (Si2^^^ and Im (Si2^^^. In this case, leptonic 
constraints, already using the present upper bound, are competitive and constrain the direction in 
which the constraint from ek is not effective (see the upper left plot). Notice that this direction is 
rotated with respect to the LL case because of the presence of LL x RR double Mis. 
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FIG. 7: Same as Fig.[T]for [^f-^^^- 

In Figure [10] we can see the bounds on Re ( ) I™ ( ) • The same bounds apply also 

to the LR case. For these Mis, the hadronic bounds come also from e' je and are quite stringent. 
However, the bounds from /i — > 67 are even more effective. Also in this case the bounds are 
independent of tan /3. 

VIII. CONCLUSIONS 

While there exists a huge literature dealing, separately, with FCNC constraints on the hadronic 
and leptonic SUSY soft breaking terms, much less attention has been devoted to the intriguing 
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FIG. 8: Same as Fig.[T]for ((5f2)LL- 

possibility that the two sectors may find correlated bounds in SUSY theories with an underlying 
grand unified symmetry. 

We have pursued such an analysis in the context of a broad class of theories which are based on 
two appealing assumptions: i) local SUSY is broken in the observable sector through gravity medi- 
ation with the corresponding soft breaking terms arising (as momentum-independent hard terms) 
at an energy scale close to the Planck mass; ii) the fundamental gauge symmetry of the theory 
includes a grand unification of quarks and leptons which is present down to the typical GUT scale 
and hence constrains the form of the supergravity Lagrangian, in particular its Kahler potential. 
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FIG. 9: Same as Fig.[T]for (<5i2)rr- 

The presence of the general conditions i) and ii) entails some correlation between hadronic and 
leptonic soft terms at the superlarge scale where they first arise. 

Obviously, the extent to which such correlation survives when performing the superlarge run- 
ning of the soft breaking terms from a scale close to the Planck mass down to the electroweak scale 
depends on the new physics present in such long interval. Here we adopted the two simplest possi- 
bilities: just a big desert or a new intermediate scale, below the GUT scale, where the right-handed 
neutrinos acquire a mass in a SUSY see-saw framework. Moreover, to make the problem treatable 
in a model independent way, we made two relevant simplifications on the possible pattern of the 
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FIG. 10: Same as Fig.[I]for {Sf2)RL- 

soft breaking terms: we considered that only one of the FC Mis in Eqs. ©-([T]) is switched on at a 
time and in the discussion of the hadronic FCNC constraints we took the gluino exchange to be the 
representative source of the SUSY FCNC contributions. On the other hand, in the leptonic sector 
where there is no analog of the gluino dominance, we performed a full computation (in the slepton 
mass eigenstate basis) pointing out the possible cancellations which may arise when the various 
SUSY contributions are taken into account. Indeed, the first part of the present work provides a 
renewed, thorough and comprehensive assessment of the bounds on the hadronic and leptonic Mis 
taking into account all the relevant pieces of information on the FCNC phenomenology we have 
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TABLE VIII: The final score of the quarks vs leptons match showing the dominance of hadronic bounds. 



been accumulating in these last years. Obviously this constitutes the basis for the subsequent work 
of correlating the hadronic and leptonic FC Mis which is the main goal of the paper. 

The extent of the impact of such correlation on the present upper bounds on the FC 6 param- 
eters has been exemplified in the study of the role of LFV processes in constraining the hadronic 
6 parameters in the down-squark sector. The relevant hadronic processes (in kaon and beauty 
physics) which are involved in bounding the (AB=LL,RR,LR,RL) already provide rather 

stringent limits (see Table IIIB on most of them. Yet, the inclusion of the correlated constraints 
arising from /j ^ + 7 proves to be extremely powerful for quite a number of such (5's. This is 
the case for ('^23)j^j^ (Fig. ED as well as for [^12)^^ and ('^12) (Fig.fTOl). 

Leptonic bounds are competitive for (5^2 ) (Fig. E]), ((5^2 ) (Fig. [9]) and (sf^) (Fig. [6]). 
Interestingly enough, this holds true even when we consider the present upper bounds on the 
relevant LFV processes and, at least in some cases, it becomes dramatic when we take into account 
the future experimental sensitivities to LFV (namely, the third column of Table HVl) . On the other 
hand, most of (sfj^ FC insertions are essentially unscathed by the inclusion of the related bounds 
from LFV. This is the case for LL, LR and RL Mis in the 13 and 23 sectors. 

The above considerations provide a precious tool in the effort to disentangle the underlying 
SUSY theory in case some SUSY particles should show up in LHC physics. It will be very difficult 
to have some "direct" signal of the presence of a grand unified supergravity (for instance, through 
the observation and study of proton decay modes). Looking at correlated SUSY contributions 
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in hadronic and leptonic FCNC processes, together with some information on the scale of the 
soft breaking sector from LHC, we could have some hints on whether there exists a grand unified 
underlying symmetry. Admittedly, even if we are particularly lucky (for instance we observe SUSY 
particles and some LFV processes), this is going to be a long term project which requires a lot of 
sweat and educated guesses. But, at least, our paper indicates possible paths to follow to achieve 
some result in the difficult task of "reconstructing" the correct fundamental SUSY theory. 

In this sense, our work is yet another relevant proof of the complementarity of flavor and LHC 
physics in shedding light on such an underlying new physics beyond the SM. 

On a more phenomenological ground, our results can find an important application in individu- 
ating for each FC 6 MI which process (either hadronic or leptonic) is more suitable to constrain it. 
If it is true that in some cases hadronic (5's find a better limit when LFV processes are taken into 
account (as we discussed above), also the reverse turns out to hold in several circumstances. For 
instance, we pointed out that there are cases when no bound emerges from LFV for some leptonic 
6. This happens for [623] (see Table IVTH) and in this case we have to make use of FCNC in B 
physics to extract a bound on such leptonic FC quantity. There is a healthy competition between 
hadronic and leptonic FCNC physics in limiting the SUSY Mis. A comprehensive score of such 
hadron versus lepton "match" is provided in Table IVIIII which shows the final ranking: quarks win 
with 21 points and leptons follow with 12 points.^ 

In conclusion, we hope that this work may display the richness which is present in flavor physics 
once we assume a grand unified supergravity framework with gravity mediated SUSY breaking. It 
could be that, at the end, flavor physics is one of the very few handles we have to understand from 
low-energy physics whether Nature has chosen to possess supersymmetry and grand unification at 
the root of its symmetries. 
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